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Abstract: This paper deals with the queuing system M/M/l with additional servers for a longer queue. 
Clearly the traffic intensity for this system will depend on the number of additional servers. The expected 
number of customers in the system, the probability of the additional of one server and the probability of 
the additional of two servers are obtained under the assumption that the number of additional servers 
depends on the number of customers in the system. The condition under which the M/M/l queuing system 
with additional servers is profitable is discussed. A MATLAB program is used to illustrate this condition 
numerically. Finally, the maximum likelihood estimators of the parameters for this queuing system are 
obtained. 
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I. Introduction 

In many queuing situations, normally the number of servers is dependent on the queue length. For 
instance, in a bank more and more windows are opened for service when the queues in front of the already open 
windows get too long. This procedure is sometimes used even in the case of airlines, buses, etc. The basic 
advantage in such a system is the inherent flexibility of service. Here we shall consider a simplified system of 
this type. 

Suppose we deal with a queuing system in which the arrivals are individually in a Poisson process with 
parameter X, and service times of customers are negative exponential with mean a" . The first customer to come 
is the first to be served. As long as there are enough customers for service. As long as the number of customers 
in the system is greater than or equal to zero and less than or equal to N, there is only one server in the system. 
As the number of customers in the system increases to more than N and is still less than or equal to 2N, an 
additional server is added. This additional server is removed when the number of customers in the system 
decreases to N or less. As soon as the number of customers in the system goes beyond 2N, the number of servers 
will be three. Similarly, the third server will be taken off when the number of customers falls to 2N or below. 
Let A, n and \i n be the arrival and service rates, respectively, when there are n customers in the system and so, for 
the described queuing model, we have 
X n =X for every n>0 , 

if 0 < u < N , 
if N < ti < 2N, 
u n = i J ^ if 2X < n . 
where, X, [i and P ~ be the arrival rate, service rate and traffic intensity respectively. 



(2n; 



II. System Equations 

Assuming the limiting distribution of the number of customers in the system to be { }. We obtain the 
following steady-state equations. 

( C+y )P lB = (P r ;n - 1) + ^(n + 1 ) ; 1 < n <N , 

(C+M )PiN = (P L (N - 1) + 2uP|(N + 1) , 

( (+2y ) P i n - (P r > - 1) + ZuP^Cii + 1) ; N<n<2N , 

( (+2iJ )P t ZN = (P (2\ - 1) + 3uPi(2 \ + 

( (+3jJ )P L n = (P r ;n - 1) + 3(^(11 + 1) , n>2N . 
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Solution of the System Equations 

This system of equations can be solved recursively, or by using the general solution for single 
Markovian queuing models which is given as follows 
For 0<n<N ; 

p n =niG = orcn - 1)1 KcGi/mc-h i) )p^o a 

= n Po : (writing M\i =P ) 

= P**Q . (2.1) 
For N<n<2N ; 

p n = ni(i = o)> - i)| K (Ci/ ma + i) )p^o a 
= [fUO = 0)^ - 1)1 IC(ti/|ti(i + 1) )3 n L (i = >J) r (n - i)| Icci/ittCi + 1) )3 ]P 0 

= (^/u) N .(mu) n - N Pq 

=2 n-N H 0 _ (2.2) 
For 2N<n ; 

p„ = n L o = orcn - 1)1 E(( L i/ ^(i + 1) )p L o a 

t 

niO = o)f(N - i)l ^(i + 1) )3 riiO = \) r (2\ - 1) [«.!/ wo + 1) )3 n^i = 2M)> - tjf K^i / + 1) )j 

= (H0 N . a/2(j) N . a/3|j) n - 2N P ° 
= 2 N 3 n " 2N pnPD . (2.3) 

00 

Now, to determine Pq we can use the normalizing condition ti=o and equations (2.1),(2.2) and 

(2.3) as follows 



N co N* 2W 

oo 



n=Q +ti=N+ 

y — —P» 



2 N 3 



N vr . i N-l , Hi1 00 

= ti=0 + ti=0 + ti=0 

i-.W 



1- P 



: ^ +1 (l-(£) ) p3 N + i 



= 1-p + 2 - p +2 w (3-p) 
1 p^ 1 p 2 ^ 1 



1 - 



(l-p>2- p) 2 N (2-p)(3-p) 



Pq" 1 = (1 - p>2 - p) 2 N (2-p)(3 -p) . ( 2.4) 

III. Profitable of using Additional Servers 

Suppose the cost to the system due to the presence of a customer is $Ci and the cost of bringing in an 
additional server is $C 2 . It is profitable to use the additional server only if its cost is less than that due to the 
waiting customers see [4]. 

Let Q be the long-run queue length of the standard queue M/M/l with the same parameters. It is evident that the 
M/M/l queuing system with additional servers is recommended only if 

C 1 [E(Q)-E(Q 1 )]>C 2 [P r (N<Q 1 <2N)+2P r (2N<Q 1 )], (3.1) 
where, 

E(Q)=Expected number of customers in the system for the queuing model M/M/l with no additional 
servers (standard). 

E(Qi)=Expected number of customers in the system for the queuing model M/M/l with additional servers. 
We have 
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E(Q) =ti=o 



2^ 



Pc 



=n=Q 

P 



= 1 - P . (3.2) 

Also, 

E(Qi)= ti=0 = ti=ij +ti=N + l +ti=2N-l 



(3.3) 

Now we have to find 

N N 



ti=0 = ti=iJ 



N 

d 



ti=0 



= (1 - p) 2 [l+N.f> N+1 -(N+1).P N ]. 



(3.4) 

Similarly, 



Si 1 - frf 

-r£i " =2 N ' °dprl 2 i 



i_£ 

n^Nil " =2 1N ' "1P[\2/ 1 2 ]= 



(2 - p) 2 P 0 [2(N+1) - 2" N+1 (2N+1)P N -~NP + 2" N+1 NP N+1 ]. (3.5) 

And 



00 

1 



n =2N-i = 2 N *P 3 - P ]= 



2 ( 3 ~ p) [3(2N+1)-2N-] . 
(3.6) Substituting by (3.4) , (3.5) and (3.6) into (3.3) it follows that 

1 

E(Q0 = .- p 0 { 0- ~ P) 2 [1+N^ -1 -(N+1)P N ] 



+ (2-p) 2 [2(N+1)- 2- N+1 (2N+l).P N -NP + 2- N+1 N* M+1 ] 

.-.2N 



+ 2 N (3 - p) 3 [3(2N+l)-2Nfl } 
= eP 0[ A/(Z T 51 E(l - p)S 2 £(2 - p)S r 2 1(3 - p)3 r 2 ) ] 

(3.7) 
where 

A= 2 N (2-P ) 2 (3-P f [1+N.^ M+1 -(N+l)- N ] 

+ 2 N (1-P ) 2 (3-P ) 2 r N [2(N+1)- 2" N+1 (2N+1)P N -NP + 2" N+1 NP N+± ] 

+ (1-P ) 2 (2-P fp 2 '' [3(2N+l)-2Nr ] . 
Also, we have to find the probability of the additional of one server, 



y /> y j— 

Pr( ;v <q 1 < 2Jv ) = " =rMi 2 "~" 
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P a 2 



7L=jV + 1 



1- £ 



(2 " Pi . 

And the probability of the additional of two servers , 



(3.8) 



Pr( 2A <Q 1)= ^jy_ 1 



.v-1 



Substituting by (3.2), (3.7), (3.8) and (3.9) into (3.1) it follows that the modified M/M/l queuing system is 
recommended only when 

C 2 E(0)-E(0 1 ) 



C 1 < r r (A < 0 1 < 2a) + 2? r (2N < 



1 - 



+■ 2/ 



1 - g)"" 0 *» 
r v \2/ , 2 £ 

' (2-p) + Z Z*0-A 



-a 



1 - p 2- (1 - p) 2 (2 - p) 2 (3 - p) 2 

(3 - p)p- v (2- v - p- v )+ 2(2 - p)p=- v / 

/2 v (2-p)(3-p) 



Substituting by (2.4), we get 



2 - 0 — p 



.v-l 



i.v-i 



Ci < [ft - p) s (2 - p) 2- v (1 - p)(2 - p)(3 - p) 2 V (1 - p)- (2 - p) 1 (3 - p) 1 ] 
2*C2-p)(3-e) 



[(3 - p)p' v (2* - p« ) + 2(2 - p)p M ] 



2 v .3-,p.(2- p-.p-'"" 1 ) 

< [ a - p) 3 



(1-p) (1 - p) 3 (2 -p)(3- p)j 



[(3 - p)p v (2 v - p*) + 2(2 - p).p i v ] 

1 2- v (2 ~ ,P)(3 - p) 2 (2 - p - p' v+1 ) - (1 - p)(2 - p)(3 - pfc™- 1 ~ A 
< 1 - ,P [ (1 - p)(2 - p)(3 - p)[(3 - p)p- v (2 V - ,p jV ) + 2(2 - p)p : - v ] ]. 
Substituting by A from equation (3.7) 



2 v (2-p)(3-.p) 3 



(2-.p-.p' v+1 )- ( 2 - P y 
(1 + .Vp''- 1 - (N + 



Ci < i - P { 

. (1 - p)(2 - p) [(3 - p) p\2N + 1] + (1 - m - p)p*2N ( 3 (2N+1) - 2N h/0] 



■2 N (1-.P ) 2 (3-P ) 2 P" V [2(N+1)- 2" N+1 (2N+1)P J> -N? +2" N+1 NP'' 1 ]} / 



■ - N+1 (2N+1)P' V 



(1 - p)(2 - p)(3 - p)[(3 - p)p v (2* - p v ) - 2(2 - p)p^ ] 
This can be simplified to 
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C t < CI - p)[(3 - pX2- v - P- V ) " 2(2 - p)p- v ] 
2" (3 - p)[2(.Y + 1) - (3.Y + 2)p + Ag s ] 

{ 0-p) 
_ 2 jV (l - pX3 - p)[2(tf + 1)- Z^^CZJV + Dp" - Np - l-x^Np* 1 * 1 ] 

_ p N [6(2N 4- 1) - 2(1LV -I- 3)p + (UN + 5)p 2 - 2Jfp a ] 

&-pi }. (3.10) 

From this inequality, depending on the decision variable, the largest value of C 2 /Ci can be determined when 
P and N are given. Table (1) has been obtained from the inequality (3.10) for the purpose of illustration. 
MATLAB program is used to perform these calculations. 

Table (1) 



fNP 


0.1 


0.2 


0.3 


0.4 


0.5 


1 


4.1765 


4.6119 


5.1762 


5.9293 


6.9778 


2 


6.5671 


7.4620 


8.5872 


10.0580 


12.0807 


3 


8.8032 


10.0205 


11.5852 


13.6639 


16.5562 


4 


11.0264 


12.5286 


14.4718 


17.0736 


20.7233 


5 


13.2486 


15.0296 


17.3346 


20.4265 


24.7769 



From this table, we can conclude that for any given value of N, the largest value of c 2 / c i decreases 
with decreasing values of P , that is with increasing values of service rates. The basic advantage in such a 
system is the inherent flexibility of service. Then there exists an optimal stationary policy (resulting in the least 
long-run expected cost) characterized by a single positive finite number N such that it is optimal to use the slow 
service rate when the number of customers in the system is less than or equal to N, faster service rate ^¥ 

when the number in the system is greater than N or less than or equal to 2N and more fast service rate when 
the number in the system is greater than 2N. 

III. Maximum Likelihood Estimators of The Parameters 

To find the likelihood function it is necessary to find the three basic components see [5]. 

1. the stationary distribution of number of units in the system with contribution ® ^ u , 

2. the interarrival times of length ti [each of which is exponential for n arrival units with contribution 



i>" 

(: = l )], 



3. the service time of durations ti for k units with contribution 
Hence, the likelihood function may be written as. 



( =- ) • 



i>- i: n 



IIS-*** 



2 jV 3 Lt_3 - v 

=:=1 :=1 ; 

(use equation (2.3)) 
JV> e\-T T) e\-uty (»(* + if} u\k - ^) k L l , 

where , 

n k 

T=i=i , t=i=a and ki is a constant . 
Thus 

ln[L{&)] = hF L 0 - (T - at -f (n -f - u~)lnu -f fnfejl. 

By partial differentiation with respect to ( and 14 , respectively, and equating the results by zeros, one can 
obtain 
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n 4- 

-CK[3ifii./saJ ^ = *rj= M =o (4.1) 

and 









6=£ = " 



5^ 



(4.2) 



where ( and QUOTE $ P- are the maximum likelihood estimators for the parameters and P , respectively. 

To obtain ( and F we need to use the following important theorem 
THEOREM (4.1), 

For the birth-death process with rates, namely 0— (i/ 1 = Ci?!* 1 , ^ — = for n positive integer and 
both S ri ana r r L are bounded functions of n [1], 
We have 

dlnP 9 



where r a = delay probability and 

E(n)= expected value of n . 
Proof: It is clear that 

Pn = W W = °?(* - i>l(i7.«ia + 1> Vn > 1 
? 0 = K[i + 5.C n = i^llLO = 0) T C" - 1)1 G ifciQ + 1) 1 a r c- 1) 

Case (i): putting tl l = C£?l z , it follows that 

P* = + 5 = E (*« n 4 (t = 0) r (?i - i>lffi^iC" + l) 3 ]3 T (-l) 

Taking the logarithm of both sides of this equation and differentiating the resulting equation partially with 
respect to C y ie ^ d 

-(Mn/> t 0)/B( = F L Q [2 L (7i = l) r ^| ln(\n - 1) U L (i = Q}\n - l^M^l) 3 ] . 

But, 

E{n) - Zi(n = l) T oo|EnF i? i = F t 0 [^(n = 0 = «) r C-i - 1) (t=/j"tC + 1) 3 ]3 . 

Therefore, it follows that 

Case (ii): putting ft = , it follows that 

? a = K[i + 5. (" = i) r «>l Ki/^ n^O = 0) r C-i - i)IC^/nO + 1)3 ]3 T (-i) . 

Taking the logarithm of both sides of this equation, it follows that 

lnP n = -^i[l + IX" = i) r co| Kl/a T n n L 0 = O) r 0i - VMii^iO + 1) 3 ] 

Therefore, 

3u = PjO [Z^n = l) f oo~ Kn/^Cn + 1) n^O = <0 T 0i - IDS'G^iP + 1) 3 ] 

Hence, 

,U 0fi =F|0 [ZiCn = 1} t ooS En/^n [JiG = 0] T (* - IjlO^iO + 1) 3 ] 

Thus, it is evident that 

r(>i) _■ ^ 

Hence, the theorem is proved. 

Using theorem (4.1), equation (4.1) can be written in the form 
-C(tffrLPj0)/3(+(T - (n+ u) = 0 , 

*~oo - a - + ^) = o , 
i 



T[(n + u) - E(n)] 

Also, by using theorem (4.1), equation (4.2) can be written in the form 
S(dlnP 0 ) 

■ — — 5t+ (k-u) = 0 
E(n) - + (fe- u) = 0 



(4.3) 
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u = — 

fc[(fc-tO +£(«)] (44) 
Where - is the maximum likelihood estimator of E(n) . 

Obtaining the value of - from each of the equations in (4.3), (4.4) and dividing the resulting values yields 
From (4.3) 

"(m) = ( n + 3"C - 

From (4.4) 

£u = (k — u) -h ZTfri) 9 
E(n) = t£ - Of - u) 9 

ftiT- (fc - «)]/C-rC + (* + m = 1 . (4.5) 

Recalling that P , the maximum likelihood estimator of the traffic intensity P , is given by P = C fp and 
using equations (4.3), (4.4) 

p' = C/^r = (t[0i + ^) - - a) + F»]) (4. 6) 

This equation can be written in the form 

t[?i + u — F(?i)] - pT[,k- — u + £"00] = 0 (4 7) 

/? 

Consider the special case M/M/l on replacing £"00 by 1 — 0 into equation (4.7) 

t n +- it - -——J - pT fc - m +- -—J = 0 
1 - p] [ 1 - p] 

t[(n -f u)(l -p)-p\- pTVji - u)(l - p)+p\ = Q 
T(k - u - l)p 2 - [t(n + u + 1) + T(k - u)]p + t(n + u) = 0 . (4.8) 
Replacing (k-u-1) and (n+u+1) by (k-u) and (n+u), respectively, in (4.8) yields. 
T(k - u)p= - [t(n + u) + r(fc - u)]p +- tfo + u) = 0 P 
TCfe - lOtfGf - 1) - + *0(p - 1) = 0 
0?" - «>p - t(n -f u)] = 0 (4 9) 

r(fc-u)p- t(r\ + u) = 0 
Then 

tfri + 



p = ■ 



(4.10) 



T(k - ui J 

is a solution for the modified equation (4.9) of degree two in P . Using (4.5) and (4.10) it follows that 
p" = t(n + E0/T(fc - * (tdT - (fc - ufift-7C + C?i - *0D , 
p~ = {^2 if(n -h - £(fc - + !£))/(- F r 2 (X L : - a) + - z0(?i + ")) . 
—7*2 pCtk - *0 + r P"C^ - u)(r\ + u) = t T 2 4- u) - t(k - u'){n 4- u] , 

t T 2 iffta + u) -h r T 2 - u) = 0 1 + «Xfe - *0 + *] . (411) 

To find P , put C = P P into equation (4. 1 1) 

t 2 S(7i + «) + r s p s ,u(*: - = (?i + - u){Tp+ t], 
(n + i0O- n)[rp+ r] 

(use equation(4.10)) 



+ + s p s (r: — u) 

t(ti -h □ ) 



(ti + u)'k — u" 



Tfk-u 



-h t 



t»6, + ») + T^(k-.) 

T * (k — u)^ 
t( ri 4- u )[ii — ii — k — n] 



t ' ti - u 
n + k 



t(ti-h u! 



k-u 



t-h T 



t' ti - u ■ 
Tk-u. 



(from equation(4.10)) 



Thus 
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A n+k 

y — 



t+Tp . (4.12) 
To find 0 put = C/f r into equation (4. 12) 
OP" = (n + k)/(Tp~+t) , 

Tp + t (413) 

Then, by using (4.10) it follows that 
n -h k k - y 

^ = TpTt = ~~T~ 

(4.14) 

f = C P (n + k))/(Tp + t) = (n + u )/T . 
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